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Abstract. In this talk I describe work on computing non-leptonic matrix elements 
consistently with both long and short distance contributions included. On the simpler 
example of the tt^-tt^ mass difference I explain in detail the matching procedure and 
^+ ■ the difference between various low-energy models. I then explain the new difficulties 

^^ \ in non-leptonic Kaon decays and how the matching here can in principle be done in 

the same way when scheme dependences are correctly accounted for. In the end I 
summarize the results J. Prades and I obtain for the A/ == 1/2 rule and Bq. 
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J": INTRODUCTION 

^ . The problem of describing non-leptonic decays is a very old one and is still not 

fully solved today. In this talk I will describe the large Nc method first suggested in 
K^ ; a series of papers by Bardeen, Buras and Gerard [1] and later extended by several 

;_, ■ other authors. I will illustrate most of the problems and solutions on the example of 



c^ 



the charged and neutral pion mass difference and afterwards show how this method 
can be extended systematically to the case of non-leptonic weak decays as well. The 
main results described there are those of [2] and also descibed in [3]. 

The subject of this meeting was hadronic physics, so why are we interested in 
these extra quantities. They provide a very strong test of our understanding of 
the strong interaction at all length scales. Our present knowledge of the strong 
interaction can be summarized as: 

Short Distance: This is the perturbative QCD domain and here QCD has had 
many successes, we count this region as understood. 
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(Very) Long Distance: This is tlie Cliiral Perturbation Theory (CHPT) regime 
[4]. Many successes again and basically understood. 

Intermediate Distance: This is the domain of models supplemented with vari- 
ous arguments, sum rules, lattice QCD results, etc. and is the most difficult. 

In the type of observables covered in this talk all three regimes are important. 
We consider processes with incoming and outgoing hadrons but with an internally 
exchanged photon or weak boson. The difficulty now resides in the fact that even if 
the external hadrons have all low momenta we need to integrate over all momenta 
of the internal 7 or W^. This means that all regimes come into play and that they 
need to be connected properly to each other. The last is known as matching. 

The main part is in Sect. I where I show how we can explain the mass difference, 
m^+ —ml.0 using this class of methods. Here we can also see how the model approach 
and the correct answer agree. Sect. II then covers the extra problems involved in 
non-leptonic weak decays and how the X-boson method of [2] can be used to solve 
those. Finally I present numerical results for this case and conclusions. 

I A SIMPLE EXAMPLE: THE 7r+-7r0 MASS 

DIFFERENCE. 

This non-leptonic matrix element has several features that make it simpler. 

1. We can neglect rriu and rud to a rather good approximation. This then allows 
current algebra to relate the electromagnetic mass difference to a vacuum to 
vacuum matrix element only [5]. This can then be related to the measured 
hadronic cross-sections in electron-positron annihilation so in this case we 
know the correct answer. 

2. There are no large masses involved so there are no large logarithms that need 
resummation. 

3. The photon itself provides for an easy identification of correct scales. 
Basically the procedure is now to evaluate 

where M stands for the meson under consideration and J^ for the electromagnetic 
current. ^ is a gauge parameter. The procedure is now as follows: 
l:We rotate the integral over photon momenta in Eq. (1) to Euclidean space. This 
has two advantages, in Euclidean space thresholds and poles are smoothed out 
making treatment of these easier and Euclidean space momenta have all components 
small if g|; is small. The latter allows for a simpler identification of long and short- 
distance than in Minkowski space. 
2: The final step is now to set 



d^QE = / qldqE dVt = q%dqE dVt+ q%dqE / dVt (2) 

JO J Jo J J a J 

' . ' ^^i! , ' 

long-distance short-distance 

and perform both integrals separately. Notice that the scale /x is just a splitting scale 
in the integral and is not directly related to any subtraction scale in the calculation 
itself. Therefore, if both the long-distance (from to ^) and the short-distance are 
calculated with high enough precision the final result should be independent of the 
value of /i. We check this by varying /i in all our calculations, i.e. we check the 
matching. 



A Short-Distance 

The short-distance contribution was first calculated in [6] using the sum rule 
by Das et al. [5]. It was later rederived using the Operator Product expansion in 
[7]. The diagrams in Fig. 1 depict the main contributions. Performing the photon 
integral leads to a set of four-quark operators that can be evaluated in leading 1/Nc 
since we can then apply factorization. The result is [6,7] 

2 2 \ Sa^a 2d2 ^asa{qqy 



[mU - <oj^^ = ^F^B^ = -^^ , (3) 

with F the pion decay constant in the chiral limit and Bq the parameter in lowest 
order CHPT describing the quark condensate. 



B Long-Distance 

In the previous subsection we could use perturbative QCD but that is not possible 
in the long distance domain. So here we have to put in the things we know and try 
various models. 
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(a) (b) (c) 

FIGURE 1. The three short-distance contributions, (a) electromagnetic quark-mass corrections, 
(b) Penguin-hke diagrams (c) Box Diagrams. The dashed hne is a gluon, the dotted hue a photon 
and the fuU hues are quarks. 
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FIGURE 2. The long-distance contributions to (m^^ 
the fuU hnes are pions. 



The dotted hne is a photon and 



CHPT: This can be done for /i rather small and leads to 
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The first term was first done in [6] and the chiral correction in [8,9]. The two 
contributing diagrams are depicted in Fig. 2. Something that is important is 
that the gauge dependence only cancels between the two diagrams in Fig. 2. 

Chiral Quark Model: This was done in [10] and gives only a marginal improve- 
ment. Note that we cannot use the usual dimensional regularization here but 
must use the cut-off in the photon propagator. There is the additional prob- 
lem that at first sight only the equivalent of the diagram of Fig 2(a) appears, 
which is a two-loop diagram, and the result is not gauge invariant. Only after 
the equivalent of (b) is added, which is a three-loop diagram, does the gauge 
dependence cancel as required [10]. 

With Vector-axial-vector Mesons: We have to include here Weinberg's con- 
straint on the couplings to obtain a unique result otherwise the result will 
be very dependent on the specific model used. E.g. a hidden gauge model 
with only vector mesons is still quadratic in /x^ but with a negative coefficient. 
Using Weinberg's constraints leads to 
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But beware of partial results. Using a linear vector representation only even 
gave a quartic dependence on fi [8]. The result in (5) for /i — > cx) is basically 
the result of [5] and was also obtained in [11]. It has also several nice features. 
Expanding in n for small /i reproduces the CHPT result with the meson dom- 
inated value for Liq. For large ^ it goes as l//i^ so it can match on very well 
to the earlier short-distance result. 

ENJL This basically coincides with the previous result and was first obtained in 

[12]. 

Extensions of the above exists for nonzero-quark masses [7,13] and references 
therein and also with more large A''^ arguments to underpin the matching [9]. 
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FIGURE 3. The long-distance contributions and short-distance contribution compared with the 
experimental values. 



C Discussion 

Numerical results are shown in Fig 3 for all cases. The experimental value and 
the one with the sub-leading in 1/A^c subtracted are shown as the horizontal lines. 
The subtracted part is the chiral logarithm contribution as estimated in [13]. Notice 
that CHPT starts deviating quickly from the VMD and ENJL results. The CHPT 
result is only reliable up to about 500 MeV. The VMD result and the ENJL result 
basically coincide here, the difference is due to the precise input values. Both these 
curves also follow essentially the correct answer as obtained from electron-positron 
annihilation and the sum-rule of [5]. 

Notice that the ENJL model has the correct matching on to the low /i CHPT 
result and is a considerable improvement over it at higher /i. Notice also the almost 
perfect agreement with the estimated part leading N^ part of the mass-difference. 

From this section we can conclude: 

1. Different Low energy models give quite different results and we have to use 
short- distance constraints and phenomenological inputs to improve the long- 
distance contribution to above the regime where CHPT is applicable. 

2. CHPT alone for the long-distance regime is as a first guestimate acceptable 
but start differing from the correct answer at a scale of about 500 MeV. 

3. Even for this low-momentum dominated observable the short-distance contri- 
butions are sizable at scales around 800 MeV. 
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FIGURE 4. The naive W^ exchange contribution to K^ -^ tt+tt'^. 

II KAON NON-LEPTONIC DECAYS 

One of the difficult unresolved problems is to understand the origin of the A/ = 
1/2 rule. The underlying process is W^+-exchange leading to an operator of the 
quark-structure {su){ud) which has both isospin 1/2 and isospin 3/2 pieces. If we 
assume the W~^ couples directly to hadrons the process K~^ -^ tt+tt^ goes simply via 
the diagrams in Fig. 4, but there are no such diagrams for K^ -^ tt^tt" because of 
charge conservation. So we would expect that r(i^+ -^ tt+tt^) 3> T{K^ -^ Tr^vr"). 
The experimental numbers are r{K~^ -^ tt+tt^) = 1.1 10~^^ MeV and r{K^ -^ 
tt^tt") = ^T{Ks — > TT^Tr") = 2.3 lO""*^^ MeV, precisely the opposite. Translated into 
isospin amplitudes for the decays, see e.g. [14] for the precise definitions, we obtain 
|^o/^2|exp ~ 22.1 .The problem is not due to chiral corrections since using the 
estimate of [14,15] we can extract them and get 

l^o/A2|chiral = 16.4 W2 . (6) 

naive 

where the last number is the one using naive VT^-exchange as depicted in Fig. 4. 
In the notation used in [2,14] we have 

Ao = Ci9G8 + G27)V6/9Foiml-ml) A2 = C lOG 27 Vg/ 9 Foiml - ml) (7) 

which after subtracting the estimated chiral corrections from experiment yields 

G8 = 6.2±0.7 Gsr = 0.48 ± 0.06 . (8) 

Both Gg and G27 are equal to one in the W^^-exchange limit, the constant C was 
chosen to have this. We thus have to explain the large deviation from 1 using the 
corrections suppressed by l/Nc. 

This is not a hopeless task as the sub-leading corrections coming from the dia- 
grams in Fig. 1 with the photon replaced by the gluon are of order 
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compared to the leading contribution and this is in fact larger than one. 

Luckily we know how to resum this type of logarithms [16]. At a high scale we 
can replace the effect of H^^-exchange by a sum of local operators by virtue of 
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FIGURE 5. The diagrams needed for the identification of the local operator Q with X-boson 
exchange in the case of only one operator and no Penguin diagrams. The wiggly line denotes 
gluons, the square the operator Q and the dashed line the X-exchange. External lines are quarks. 

the operator product expansion. We can then use the whole renormahzation group 
machinery to run this sum over local four-quark operators down to a low scale /i/j. 
This is explained in great detail in [16]. The end result is 

Hw= T. C^{fiR)Q^ (10) 

j=l,10 

with a series of known coefficients Ci{fifi), the Wilson coefficients. The final answer 
is then the matrix element of this sum over four-quark operators, {7nr\Hy\r\K) . 
But here we have two problems: 

1. In the previous section the short and long-distance contributions were sepa- 
rated via the photon momentum. Here we have to link this somehow to the 
scale ^R appearing in the weak Hamiltonian Hy\r. 

2. To next-to- leading order in the renormahzation group the coefficients Ci{fiR) 
also depend rather strongly on the precise definition of the local four quark 
operators Qi in QCD perturbation theory. 

In [18] we showed how the method of [1] supplemented with the correct mo- 
mentum routing [6,17] solved problem 1. In [2] and in the various already pub- 
lished talks [3] we showed how a careful identification across the long-short-distance 
boundary is also possible in this case. The basic idea is to go at the scale fiR back 
from the local four-quark operators to the exchange of a series of X-bosons. These 
X-bosons can then be treated in exactly the same way as we did the photon in the 
previous section, thus allowing a correct calculation at all length scales. 

So we replace, using a single operator as an example 

CiQi = {sLli,dL){uLlt,UL) ^^ X^ [giisLl'^dL} + g2{uLl^UL)] ■ (H) 

Using the tree level diagrams of Fig. 5 this gives Ci = gig2/ M^- If we now include 
the one loop diagrams we obtain instead: 



—2- 1 + (Xs{li)ax + asiHRJh log —^ 



Ci (1 -h as(/^ij)n) = -irjT 1 + "5(/^)ai + asil^RJh log -^- . (12) 



On the l.h.s. the scheme dependence disappears but there is a dependence in ri on 
the choice of external states. The exact same dependence in ai cancels this. 

We now split the integral over the X-boson momentum as in the previous section 

/ dpx =^ / dpx + / dpx (13) 

Jo Jo Jfi 

In the final answer all Mx dependence drops out, the logarithm proportional to 
bi shows up in precisely the same way in the evaluation of the short distance 
part of (13) which is proportional to gig2/Mx {as{fi)a2 + as{fJ')bi^og{Mx/fi'^)} 
The coefficients ri, ai and 02 give the corrections to the naive l/A'c- method. 

We now use the X-boson method described above and put ji = hr. The low 
energy part can be calculated using CHPT, this is the approach used originally 
by [1] and presently pursued by [19] without including the corrections due to the 
change in scheme when going to the long-distance part. Their results coincide with 
ourss when we restrict our results to their approximations. We obtain [2] 

Bk = 0.6—0.8 B^ = 0.25—0.40 Gg = 4.3—0.75 
G27 = 0.25—0.40 Gg = 0.8—1.1 Be{fi = 0.8 GeV) ^ 2.2 (14) 

Bk is the bag-parameter relevant for K^ — K^ mixing at the physical quark masses 
and B^ the same in the chiral limit. The quark mass corrections are quite sizable. 
The results for Gg and G27 are obtained without any free input and agree within 
the uncertainties of the method with the experimental values. We conclude that 
we basically now have a first principle understanding of the A J =1/2 rule. We 
discuss the various contributions below. Gg is the coefficient of the weak mass term, 
it contributes at leading order to processes like K^ -^ 77 [14] and is often forgotten 
in those analyses. Finally B^ is much larger than used in all the analysises of the 
recent experimental results for e'/e [20] which is very encouraging. 

The final result for Gg is depicted in Fig. 6. We have shown the one-loop 
result (1-loop), the two-loop result with NDR Wilson coefficients (2-loop) and the 
two-loop results with correction for the long-distance scheme added (SI) using our 
results for the long distance part. We also showed what naive factorization with 
the SI Wilson coefficients would give and what the method of [19] would give in 
the chiral limit with the same Wilson coefficients (SI quad). 

If we look at the various contributions to Gg we see in Fig. 7 that the contribution 
of Qi and Q2 are both large and fairly constant while Qq contributes 20% or 
less. If we look inside the calculation we see that the difference with the G27 
evolution is mainly given by the long-distance Penguin-like contributions to Q2- 
The behaviour of -Be is more difficult, it is ill-defined in the chiral limit in the 
factorizable approximation [2] and we can thus only define it with respect to the full 
large N,. limit. Calculating it in CHPT only then gives fairly low values as is visible 
in the second line of Table 1. Adding higher order corrections we immediately 
obtain a strongly enhanced value as is obvious from the second line in Table 1. 
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III CONCLUSIONS 



a The X-boson method in combination with large Nc arguments allows to 
correctly identify quantities across theory boundaries assuming we can identify 
currents across the boundary. 

b The mass difference 'm'i-+ — m^o is well described by these methods with a 
surprisingly large short- distance contribution. 

C The A J = 1/2 rule is now quantitatively understood to about 30% with NO 
free input. This calculation passes all requirements usually asked of in this context 
but there are many technical subtleties. 

d -Be ~ 2.2 is good news for those trying to explain the observed values of 
e'/e within the standard model. 

e This program has been quite successful but we need new ideas to calculate 
more complex processes. 
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FIGURE 6. The octet coefficient Gs as a function of /i using the ENJL model and the one-loop 
Wilson coefficients, the 2-loop ones and those including the ri (SI). In the latter case also the 
factorization (SI fact) and the approach of [19] (SI quad) are shown. 
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FIGURE 7. The composition of Gs as a function of /i. Shown are Q2, Qi + Q2, Qi 
and aU 6 Qi. The coefficients ri are included in the Wilson coefficients. 
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